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Abstract. Mnlti-wavelength observations made in the last decade sug¬ 
gest that the universe underwent an intense phase of star formation in 
the past (z > 1). This intensive activity is commonly attributed to a 
higher galaxy merger rate when the universe was a fraction of its present 
age. We examine the role of the powerful radio sources whose comov¬ 
ing density is known to be a few orders of magnitude higher at z ~ 2, 
the ‘quasar era’. Taking into account recent models for the temporal 
evolution of the size and luminosity of a powerful double radio source, 
as well as ACDM simulations of the cosmic web of baryonic material 
at different redshifts, we argue that during the quasar era a high frac¬ 
tion of the volume of the web was occupied by the lobes of double radio 
sources. Widespread compression of protostellar clouds, triggered by the 
high pressure of the synchrotron plasma of the radio lobes, can thus be 
expected to have played a significant role in the global star formation 
history of the universe. These lobes can also yield a rather high level of 
magnetization of the intergalactic medium at these early cosmic epochs. 


1. Introduction 

Flux density-limited samples of radio sources show that the comoving space 
density of powerful radio galaxies (RGs) declined between 100 and 1000 times 
from redshifts of 2-3 to the present epoch (e.g. Willott et al. 2001). The star 
formation rate also peaked in roughly the same epoch (e.g. Archibald et al. 2001). 
Because of severe adiabatic and inverse Compton losses, most old and large radio 
galaxies are very difficult to detect in surveys and only the youngest can be 
seen at high redshifts (e.g. Blundell et al. 1999, hereafter BRW). Cosmological 
simulations have indicated that most of the matter that will form galaxies by 
the current epoch was in the form of filaments that filled only a small portion 
of the universe at those redshifts (e.g. Cen &: Ostriker 1999). Together, these 
facts lead us to conclude that the formation of many of those galaxies may have 
been triggered by overpressured radio lobes, which probably filled a substantial 
portion of those filaments then. A preliminary discussion of this work is in 
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Gopal-Krishna & Wiita (2001, hereafter GKW) and more extensive calculations 
are underway (Osterman, Wiita, Gopal-Krishna & Kulkarni, in preparation). 


2. Radio galaxy visibility 


All modern models of RG evolution (Kaiser et al. 1997; BRW; Manolakou & 
Kirk 2002) agree that radio flux declines dramatically with increasing source size 
(adiabatic losses) and with z (inverse Compton losses off the cosmic background 
radiation). Theoretical distributions of RG powers, sizes, redshifts and spectral 
indices can be nicely matched by models that require most RGs to remain active 
for T ~ 5 X 10®y and to have a distribution of jet powers (Qo) that goes as ~ 
Qq^'^ (BRW). X-ray observations indicate that the density of the matter through 
which the jets propagate declines with distance roughly as n{r) = nQ{r, 
with no = 1.0 X 10^m“^, oq = 10 kpc, and f3 = 1.5. This leads to the linear size 
of the RG being given by 


D{t) = 3.6ao ^ 


t^Qo \i/(5-/3) 


agrupno 


) 


( 1 ) 


Using these models, we find that in most RGs, particularly those at z > 2, 
the central engines remain active for much longer times than those galaxies are 
detected in flux limited surveys, and therefore they should grow to very large 
linear sizes (typically D{T) > 1 Mpc), although detecting them would require 
extremely sensitive radio surveys with redshifts. From BRW we find that the 
visibility time, r oc , and to properly estimate the actual number of RGs 
from those detectable in flux-limited surveys, one must multiply by a correction 
factor (T/r) of roughly 50 for powerful RGs during the quasar era (GKW). 


3. Radio luminosity function — RLF 

Most RLF studies are plagued by uncertainties resulting from incomplete knowl¬ 
edge of the redshifts of the radio sources, but results based upon the 3GRR, 6 GE 
and 7CRS surveys of different flux limits have the advantage of having 96% of 
their redshifts known (Willott et al. 2001). In addition, their selection at low 
frequencies minimizes the bias due to relativistic beaming. 

The powerful FR II sources are nearly 3 dex above the local RLF by z ~ 2, 
and their RLF varies little out to the beginning of the quasar era at z ~ 3, while 
it appears to decline at higher z (Willott et al. 2001). Furthermore, the RLF for 
those redshifts is nearly flat for over a decade in radio power above P 151 > 10 ^®'® 
W Hz“^ sr“^, which is where the FR II sources are most numerous. 

Combining these results with the correction factor discussed in §2 we find 
that at 2 : = 2.5 the actual proper density of powerful radio sources born in an 
interval T is p ~ 4 x 10“®(1 -|- 2 ;)^r 5 Mpc“^(AlogPi 5 i)“^, where = T/{5 x 
10®y). We then integrate over the roughly 1.25 dex of the peak of the RLF. 
Finally, we take into account the fact that several generations of RGs will be 
born and will die within the ~ 2 Gyr duration of the quasar era. This leads 
us to the total proper density, <h, of intrinsically powerful radio sources: ^ = 
7.7 X 10“^ Mpc“^, which is independent of the assumed value of T. 
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4. The relevant universe 

Numerical models of the evolution of ACDM universes indicate that at z ~ 0, 
roughly 70% of baryons are in a cosmic web of filaments of warm-hot gas and 
embedded galaxies and clusters that together occupy only about 10% of the 
volume of the universe (e.g. Cen & Ostriker 1999; Dave et al. 2001). But at 
z ~ 2.5 the growing network of filaments comprised only about 20% of the 
baryonic mass, and a quite small fraction, rj ~ 0.03, of the total volume. 

The massive galaxies that harbor supermassive black holes large enough to 
form RGs at early times would have typically formed in the densest portions of 
those filaments. The lobes ejected from them would mostly remain within the 
filaments, and since it is in this relatively small, ‘relevant universe’, that new 
galaxies formed out of denser gas clumps, we only need to be concerned with 
what fraction of this relevant universe the lobes permeated. We find that the 
mean volume of a radio source is (V(T)) ~ Mpc^, and thus, the volume 

fraction of the relevant universe which radio lobes born during the quasar era 
cumulatively swept through is: 

( = <p (V(5 X 10®yr)) (0.03/rj)(5/RTf ^ 0.5, (2) 

where Rt ~ 5 is the typical width to length ratio of an RG. The energy density 
injected by the lobes into the filaments is ii ~ 2 x 10“^® J m“^ for those same 
canonical parameters. Because (V(T)} is a sensitive function of T, if the typical 
RG lifetime is < 10® yr then ( < 0.01 and u < 9 x 10“^® J m“®(GKW). 

5. Radio lobe triggered star formation 

The discovery of the alignment effect between extended optical emission lines 
and radio lobe directions (e.g. McCarthy et al. 1987; Chambers et al. 1987) 
quickly led to calculations (e.g. Begelman & Cioffi 1989; Rees 1989) that indi¬ 
cated that star formation could be triggered by these expanding overpressured 
lobes. Recent hydro dynamical simulations including cooling (Mellema et al. 
2002) confirm that this is likely to occur through cloud fragmentation, cool¬ 
ing and compression. HST observations of high-z RGs and associated optical 
emission (e.g. Best et al. 1996) support this scenario. 

We further estimate (GKW) that these powerful RGs create lobes which 
typically remain rapidly expanding, with overpressures of factors exceeding 100 
(or Mach numbers above 10) out to distances of well over 1 Mpc. Supersonic 
expansion into a two-phase circumgalactic medium will compress many of the 
cooler gas clumps, rapidly reducing the Jeans mass by factors of 10-100 and 
thereby triggering starbursts (Rees 1989; Mellema et al. 2002). 

6. Conclusions 

Although the local universe is very sparsely populated by powerful radio sources, 
several large factors work in the same direction to make them remarkably impor¬ 
tant for galaxy formation during the quasar era. First, their comoving density 
was roughly 1000 times higher at ~ 2 < 2 : <~ 3. Second, only a small fraction 
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(roughly two percent) of the powerful sources present during that period are 
detected in the surveys used to produce the RLFs, because of severe inverse 
Compton and adiabatic losses; these unseen radio lobes fill very large volumes. 
Third, the fraction of the volume of the universe occupied by the material dur¬ 
ing the quasar epoch that would finally condense into clusters of galaxies was 
only a few percent, so these lobes only had to permeate this ‘relevant universe’ 
rather than the entire universe. In that the best models of RG evolution indicate 
that the lobes are overpressured and supersonically expanding into the relevant 
universe, the scenario that many massive starbursts, and even many galaxies, 
are formed in this fashion, may turn out to be realistic. 

An exciting implication of this scenario is that RGs inject a substantial 
amount of energy in the form of magnetic fields into the cosmic web of filamen¬ 
tary IGM at 2 ~ 2. Equipartition field strengths of 10“^G can arise (GKW), and 
are in accord with recent rotation measure based estimates for fields within the 
denser parts of the IGM (Ryu et al. 1988). It is well worth noting that two very 
recent independently advanced arguments, starting from essentially orthogonal 
evidence based on overall energetics (Kronberg et al. 2001; Furlanetto & Loeb 
2001) came to very similar conclusions about the importance of RGs for the 
injection of magnetic energy and for galaxy and structure formation. It is clear 
that our ongoing investigations of aspects of this scenario will be worthwhile. 
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